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ABSTRACT 

We have used the University of Tasmania's 30-m radio telescope at Ceduna in South 
Australia to regularly monitor the flux density of a number of southern blazars. We 
report the detection of an annual cycle in the variability timescale of the centimetre 
radio emission of PKS B1622-253. Observations of PKSB1519-273 over a period of 
nearly two years confirm the presence of an annual cycle in the variability timescale 
in that source. These observations prove that interstellar scintillation is the principal 
cause of inter-day variability at radio wavelengths in these sources. The best-fit annual 
cycle model for both sources implies a high degree of anisotropy in the scattering screen 
and that it has a large velocity offset with respect to the Local Standard of Rest. This 
is consistent with a greater screen distance for these "slow" IDV sources than for rapid 
scintillators such as PKSB0405-385 or J1819+3845. 

Key words: Scattering - ISM: general - galaxies: active - galaxies: nuclei - radio 
continuum: general 



1 INTRODUCTION 

Active galactic nuclei (AGN) exhibit variability in their 
emission over all wavelength ranges and on a variety of dif- 
ferent timescales. At radio wavelengths the dominant emis- 
sion from AGN is from relativistic jets of material ejected 
close to our line of sight. Inverse Compton scattering lim- 
its the br ightness temperature of the radio emission to be 
<10 12 K (|Kellerman fe Paulinv-Totbl fl969). which in turn 
sets a lower limit on the angular size of the radio emission 
region. For most AGN the linear size of the radio emission 
region inferred from these arguments is of the order of light- 
weeks to light-months. Hence causality arguments require 
that any variability in the radio emission from the source 
must have characteristic timescales of weeks-months or 
longer. However, a significant number of AGN have been ob- 
served to vary at centime tre radio wavelengths on timescales 
of a few days or less (e.g. Hecschcn 1984; Quirrcnbac h et ail 
1 19921 : iKedziora-Chudczer et al.ll200ll ; lLovell et all 120031 ) . If 
this variability cannot be intrinsic to the source, it must 
be a propagation effect imposed by the medium between 
the source a nd o bserver - interstellar scintillation (ISS). 
I Lovell et al.l (|2008l ) have demonstrated that even when fac- 
tors such as Doppler boosting are taken into account, 
any AGN small enough to show intrinsic variability on a 
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timescale of days must be of sufficiently small angular size 
to also exhibit interstellar scintillation. 

The ionized interstellar medium (ISM) scatters the ra- 
diation from the source, and if the angular size of the source 
is comparable to the angular size of the density fluctuations 
in the ISM, interference effects produce variations in the 
intensity of the radio emission at the location of the ob- 
server. Interstellar scintillation is typically modele d as oc- 
curring in a single, t hin screen of turbulent plasma l|Rickettl 
ll99d ; lNarav"anlll992ft . The relative transverse motion of the 
screen and the Earth means that that the modulated inten- 
sity pattern observed at a particular location changes with 
time, causing the observed variability in the centimetre ra- 
dio emission. There are two definitive tests for ISS as the 
cause of radio variability: an annual cycle in the variabil- 
ity timescale and a time delay in the intensity variations 
observed at two widely separated telescopes. 

The variability timescale is set by the relative trans- 
verse motion of the Earth and the scattering screen. As the 
Earth orbits the Sun the projection of its velocity vector 
onto the scattering screen will change throughout the year 
and the change in the relative transverse velocity is around 
30 km s . Provided that the scattering screen does not 
have a large intrinsic transverse velocity this will result 
in significant changes in the variability timescale over the 
course of a year, which will repeat annually: an annual cycle. 
Annual cycles have been measured i n a h an dful of AGN: 
B0917+624 jjauncev fc MacquartJ l200ll : iRickett et ail 
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l200ll ), J1819+3 845 dDennett-Thorpe fc de Bruvnl l2003h , 
PKSB1257-326 jBignall et al.ll2003l ) and PKSB1519-27 3 
ijjauncev et all 120031 ) J1128+5925 ijcabanvi et all l2007h . 
The spatial scale of the intensity pattern is typically much 
greater than the size of the Earth and so the intensity 
variations observed at widely separated telescopes are 
highly correlated. However, if the separation of the tele- 
scopes is such that a change in intensity can be detected 
on timescales less than the time the changing intensity 
pattern takes to travel between the telescopes (which is a 
function of the projection of the relative velocity vector 
onto the baseline vector and the baseline length), then 
a time-delay in the intensity patterns can be observed. 
Time-delays have also been measured in a handful of 
AGN: PKSB0405-385 (fjauncev et all l2000l ), J1819+3845 
|Dennett-Thorpe fc de Bruvnl 12002! ') and PKSB1257-326 
pignall et all l2006h . Measurements of ISS induced vari- 
ability in radio sources can be used to study both the 
source and the interstellar m edium, through the te c hniqu e 
of Earth-Orbit Synthesis (|Macquart fc Jauncevl l2002h . 
From the observations made to date, the distinguishing 
characteristic of the most rapidly varying sources appears 
to be that the scattering screen is located close (within a 
few 10s of pc) to the Earth. 

Most of the heavily studied scintillating AGN have char- 
acteristic timescales for the intensity variations of the order 
of a few hours. This means that in most cases the timescale 
can be measured to reasonable accuracy in observations of 
duration 12 hours or more. The largest survey to date for 
ISS in AGN is the Micro- Ar csecond Scintillati on-Induced 
Variability (MASIV) Survey (|Lovell et al.l 120031 ) where the 
VLA was used to monitor 443 compact flat spectrum ra- 
dio sources in four epochs of 3 to 4 days each at 5 GHz. 
Variability was detected in 56% of sources but only 12% 
showed ISS at every epoch and the great maj ority of vari- 
abilit y timescales were of order 1 day or longer (|Lovell et al.l 
2008). The more typical scintillating AGN therefore shows 
variability on longer timescales and does so intermittently. 
For a source whose angular size 8 SIC exceeds the Fresnel angle 
(i.e. 9 Brc > y/X/2wD, where A is the wavelength of the obser- 
vations and D the distance from the Earth to the scattering 
screen), the size of the scintillation pattern is governed by 
the source size, and has a physical scale D9 src at the Earth, 
so that the timescale on which the pattern advects across 
an observer is T char = D0 src /Krans (where Krans is the rel- 
ative transverse velocity of the screen and observer). Hence 
the observed scintillation timescale increases l inearly with 
the eff ective distance to the scattering material. lLovell et alj 
(2008) have also shown that the scintillation timescales ob- 
served for the majority of sources in the MASIV sample is 
inconsistent with the variability being produced in a screen 
moving with the local standard of rest. So the observation 
of long timescales combined with evidence that the major- 
ity of screens that cause ISS are not moving with the local 
standard of rest, together suggest that scattering is likely oc- 
curring at distances from Earth of many hundreds of parsecs 
in most cases. 

The long characteristic timescales prevalent in the ma- 
jority of IDV sources presents some problems for their study. 
It makes measurement of a time delay difficult or impos- 
sible, since the time taken for an inflection in the inten- 
sity is large compared to the time taken for the interfer- 



ence pattern to travel the distance between the telescopes. 
While, characterising the timescale of the variations, and 
how that changes throughout the year requires near contin- 
uous monitoring. The high subscription rates for large in- 
terferometer arrays where most previous IDV studies have 
been conducted (VLA, ATCA, WRST), are such that it is 
not possible to get sufficient time on these f acilities to under- 
take such a project. The COSMIC project l|McCulloch et all 
has been designed specifically to fill this niche, using 
a 30-m radio telescope to undertaken high-cadence moni- 
toring of a small number of strong, longer timescale IDV 
sources. Here we report the results of the first two years of 
monitoring for two of the COSMIC sources PKS B1519-273 
and PKSB1622-253. 



2 OBSERVATIONS & DATA REDUCTION 

We have used the University of Tasmania 30-m radiote- 
lescope near Ceduna, South Australia to monitor the in- 
tensity of PKSB1519-273 and PKSB1622-253 at 6.7 GHz 
for a period of approximately two years. The observations 
were made as part of the Continuous Single-dish Monitoring 
of Intraday variability at Ceduna (COSMIC) project. The 
motivation, background and obs ervation techniques for t he 
project are described in detail in lMcCulloch et al.1 (|2005T ). 

The intensity of PKSB1519-273 and PKSB1622-253 
were measured using a series of four scans, two driving the 
telescope in right-ascension and two driving the telescope 
in declination. Observations of these two sources were inter- 
leaved with observations of the ca librator 3C227 and an other 
variable source AOB0235+164 (|Senkbeil et al.1 120081 ). The 
measured intensities were corrected for antenna pointing 
and gain-elevation eff ects using the dat a redu ction methods 
described in detail in iMcCulloch et al.1 (|2005l ). The sources 
observed in the COSMIC project are split into two groups, 
those with declinations north of the zenith (approximately - 
32°) at Ceduna and those which pass south of the zenith. We 
observe one group of sources for periods of around 10-15 days 
and then swap to the other group. Figure [T] shows the inten- 
sity of the two COSMIC calibrators (3C227 for the northern 
sources and PKSB1934-638 for the southern sources) over 
the period March 2003 - January 2005. Prior to each scan a 
noise diode was switched and compared to a 1 dB step. Com- 
parison of the height of the noise diode to a source of known 
(or assumed) flux density then allows conversion of the in- 
tensity measurements from VFC (voltage to frequency con- 
verter) counts into Jy. Figure[T]shows the results of applying 
a single conversion factor, determined from the first month 
of observations of PKS B1934-638 to the entire dataset. Any 
changes in the amplitude of the noise diode will then cause 
apparent variations in the calibrator sources and this is the 
cause of the small long-timescale variations observed in Fig- 
ure!]] The mean flux densities of calibrators PKS B1934-638 
and 3C227 measured from this data are 3.952 ± 0.093 Jy and 
1.925 ± 0.048 Jy respectively, agreeing wit h the known cal - 
ibrator flux densities at 6 7 GH z of 3.92 Jy (jRevnoldsll 19941 ) 
and 1.90 Jy l|Baars et al.ll 19771 ) respectively. Figure Q] shows 
a polynomial fit to the daily mean intensity of PKS B1934- 
638 and 3C227. These fits are used to determine the correc- 
tion factor necessary to scale the mean flux density (for each 
observing block) to the assumed flux density for each cali- 
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brator. This correction factor never exceeds 3.5 percent, and 
the measured RMS values for PKS B 1934-638 and 3C227 af- 
ter its application are 0.080 Jy and 0.042 Jy respectively. 

The changes in the amplitude of the noise diode ap- 
pear to be correlated with the ambient temperature. Such 
changes occur on a wide range of timescales, not just the 
longer timescales apparent in Figure [T] and we do observe 
smaller amplitude (~ ±0.15 Jy) variations in the calibrators 
on diurnal timescales (Figure [2}. A number of approaches 
have been tested to try and correct for this effect, however, 
because the noise diode is housed within the telescope struc- 
ture there is a variable time lag between a change in the 
ambient air temperature (which is measured), and the sub- 
sequent change in the temperature of the air around the 
noise diode. The temperature of the noise diode is now 
measured, but was not at the time of the observations re- 
ported here. During the summer months air conditioning 
is used to reduce excessive ambient air temperature vari- 
ations around the receiver and back-end systems, but the 
need for proper temperature stabilization of these areas is 
now recognised. The systematic variations in the calibrator 
intensity on timescales longer than a day are determined 
by normalizing the flux density, by dividing by its mean 
over the observing period and then fitting a polynomial to 
the data. The order of the polynomial is approximately the 
same as the number of days in the observing period, so this 
approach is only able to correct for systematic variations 
on timescales of the order of a day or longer. To remove 
the shorter timescale (higher frequency) systematic varia- 
tions a polynomial low-pass filtering process is applied to 
the target-source data in the time domain. An example of 
this process is shown in Figure |3]for 10 days of PKSB1622- 
253 data. An in depth description of the method, including 
the steps which are undertaken to make sure that the poly- 
nomial f it beha v es at the two ends of the time range, is 
given in ICarterl (|2008T ). This two-pass approach to remov- 
ing/reducing the systematic variations in the data has been 
extensively tested using both real and simulated data and 
demonstrated to be both robust and effective. Approaches 
which use the measured ambient temperature and a model of 
the changes in the noise diode power with temperature have 
also been investigated and, while they yield similar results, 
they are significantly more difficult to implement. Ideally the 
calibrator measurements would be used directly to measure 
the variations in the target sources, however, the different 
right ascension of the calibrator 3C227 (9 hours right ascen- 
sion) and PKSB1519-273, PKSB1622-253 means there is 
very little overlap and this approach cannot be used. 

We have compared the observed scatter in the calibra- 
tor sources (after correction for systematic effects), with ex- 
pectations from the radiometer equation and find that they 
are significantly higher. Such effects are common in contin- 
uum measurements made with single-dish radio telescopes 
and are thought to be due to small and rapid variations 
in the electronic gain of the system. There is no easy way 
to fix this problem, indeed it is one of the fundamental 
reasons that interferometers are superior for radio inten- 
sity measurements - the electronic gain variations are in- 
dependent for a pair of antennas, hence they do not corre- 
late and interferometers are typically able to achieve close 
to their theoretical sensitivity limits. We observe deviations 
from thermal noise behaviour on timescales longer than 0.1 
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Figure 2. The intensity variations observed in PKS B1622— 253 
over the period 2004 DOY 88 - 98. Systematic variations of ampli- 
tude ± 0.15 Jy are present on timescales less than one day, how- 
ever, the ISS induced variations are larger amplitude and longer 
timescale. 



seconds and the measured uncertainty in our flux density 
measurements on a timescale of 5 seconds (approximately 
the time it takes to scan across the source) is approximately 
2.5 times the theoret i cal ra diometer estimates. As outlined 
bv iMcCulloch et al l (gooi) the accuracy of the flux den- 
sity measurements is a function of the intensity, but for 
PKSB1519-273 and PKSB1622-253 over the period 2003 
- 2005 it was typically between 1 and 2 percent of the mean 
flux density (~ 30 mjy). 

The longer timescale of the variations in PKSB1519- 
273 and PKS B 1622-253 and the use of a single-dish ra- 
dio telescope rather than an interferometer has required the 
utilisation of different methods of measuring the variability 
timescale than those most commonly used for IDV stud- 
ies. Most determinations of the characteristic timescales for 
IDV sources calculate th e discrete auto-correlation function 
l|Edelson fc Krolikl li"988h or DACF, and then measure the 
width at which the correlation drops to a half (to. 5), or 1/e 
(Ti/ e ). This approach does not work very well for the COS- 
MIC data due to the regular sampling (the source rise/set 
time at the telescope imposes a regular gap in observations 
each sidereal day). When the characteristic timescale of the 
variations is close to a multiple of a day the artifacts in the 
DACF make it difficult to accurately estimate. The Fourier 
Transform of the DACF is known as the power spectral den- 
sity (or PSD) and this is much less affected by the artifacts 
than the DACF. In most cases the PSD shows a single dom- 
inant peak (see for example Figure [4} and the corresponding 
timescale is the period we have taken as characteristic of the 
variability in this study. 

Interstellar scintillation is an inherently stochastic pro- 
cess and so for some observing periods there is no well- 
defined quasi-period for the timescale of the variations. As 
a cross-check on the characteristic timescale obtained from 
the PSD we used two secondary approaches to estimate the 
variability timescale. The first cross-check was to systemat- 
ically fold the data at a range of different periods and then 
fit a simple sinusoid to the folded data. The period for which 
the amplitude of the fitted sinusoid is a maximum and the 
residuals to the fit are a minimum is then taken as an esti- 
mate of the quasi-period. The second cross-check (which we 
refer to as scintle-counting), was to locate by eye the times 
of the turning points in the intensity variation time series 
and then estimate the quasi-period by dividing the duration 
between the first and last turning points by the number of 
turning points. Only observing sessions for which the PSD, 
data-folding and scintle-counting all returned consistent es- 
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Figure 1. Intensity of the COSMIC calibrators over the period 2003 - 2005. The extremes in the scatter of individual flux density 
measurements can be seen from the spread in the measured intensity, as can the slow changes in the mean intensity, believed to be due 
to thermally induced variations in the output power of the calibration noise diode. The red lines mark the start of the 2004 and 2005 
calendar years, and a number of periods of poor data quality are labelled. 
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Figure 3. An example of the polynomial filtering used to remove systematic variations from PKSB1519-273 and PKSB1622-235 data. 
The top plot shows 10 days of raw data for PKS B1622-253, the blue line is the fitted low-pass polynomial filter (14th order polynomial 
in this case) and the red line a first order polynomial, which is fitted to the data to estimate intrinsic variations in the mean flux density 
of the source over the observing period. The second plot shows the mean-subtracted calibrator (3C227) data for the same period, to 
which the same order polynomial has been fitted as for PKS B1622-253, this is used to estimate the systematic changes on timescales 
longer than a day. The third plot shows the post-filtered, mean-subtracted flux density measurements for PKS B1622-253. It is this data 
which is used for the time-series analysis. The bottom plot shows the difference between the raw and filtered data, demonstrating the 
presence of systematic variations of amplitude ~ ±0.15 Jy on timescales shorter than a day. 



Annual Cycles m 1519-273 and 1622-253 5 





Figure 4. The discrete auto-correlation function (DACF) and 
power spectral density (PSD) for PKSB1519-273 for the session 
2003 DOY 204 - 215. The PSD shows a clear peak at a timescale 
of 0.38 days -1 , corresponding to a quasi-period in the variations 
of approximately 2.6 days. 
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Figure 6. The intensity of PKS B1622-253 versus time for a typ- 
ical ten-day period. In this case 2003 DOY 184 - 194. 



timates for the quasi-period were included when determining 
if the characteristic timescale of the variations exhibits an 
annual cycle in these sources. 
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Figure 7. There is no correlation between the modulation index 
and mean flux density in the PKSB1622-253 data. 



in the measured flux density and S is the mean flux den- 
sity over the observing period in question). This represents 
a lower-limit to the modulation index for the scintillation 
as some of the emission is expected to be from source com- 
ponents which are too large to exhibit ISS. We investigated 
how the modulation index changed as the mean flux density 
of the sources changed and found a loose correlation between 
the two quantities for PKS B1519-273, which suggests that 
a significant fraction of the change in the mean flux den- 
sity was due to evolution of the compact, scintillating core. 
In contrast for PKS B 1622-253 there is no good correlation 
(Figure [7]), suggesting that the flux density changes include 
significant contributions from non-scintillating components 
of the source. 



3 RESULTS 

Figure [5] shows the Julian-day Intensity Monitoring 
plots (JIM-plots) for PKSB1519-273 and PKSB1622-253 
over the period 2003—2005. The mean flux density of 
PKS B1519-273 at 6.7 GHz is seen to vary between about 
2.5 Jy in mid-2003 and again in mid-2004, down to approxi- 
mately 1 Jy by early 2005. PKS B1622-253 shows even larger 
variations in mean flux density, peaking in excess of 4 Jy in 
early 2004 down to around 2 Jy later the same year. The 
variability observed on timescales of months to years is be- 
lieved to be intrinsic to the source due to evolution in the 
parsec scale jets. Variations of this type are commonly ob- 
served in the centimetre wavelength emission of AGN. It is 
also noticeable that the scatter in the measured intensity for 
PKSB1519-273 and PKSB1622-253 is significantly greater 
than observed for 3C227. This is due to the ISS-induced 
short timescale variations, which are compressed by the large 
time range displayed in Figure [5] Figure [6] shows a typical 
ten-day period for PKS B1622-253, where the character of 
the shorter timescale variability can be clearly seen. 

The modulation indices (fj, = §) for PKSB1519-273 
and PKS 1622-253 over the two year observing period were 
in the range 0.02 - 0.07 on timescales of 10 days, after cor- 
rection for systematic effects (<r is the standard deviation 



3.1 Annual cycles 

Figure [8] shows the changes in the measured characteris- 
tic timescale as a function of the day of the year for both 
PKSB1519-273 and PKSB1622-253. The blue circles are 
timescales measured during the first year of observations 
(2003), while the red crosses are from 2004. Figure [8] clearly 
shows that the timescale of the variability in both of these 
sources changes with time in a manner which repeats on an 
annual cycle. Such behaviour is an unambiguous indication 
that the short timescale variability observed in these sources 
is caused by interstellar scintillation. The annual cycles dis- 
played by the two sources are similar in their overall charac- 
ter as they are relatively close together on the sky and the 
ecliptic coordinates of the source play a major role in deter- 
mining the times during the year when the variability slows 
and speeds up. Although the two annual cycles are broadly 
similar, the actual timescales measured in each epoch typ- 
ically differ by a factor of around two, with the variations 
in PKS B1519-273 generally being faster. The observations 
used to measure the timescale in these two sources were 
interleaved and so the detection of different values demon- 
strates that they are not due to residual systematic effects. 
Timescale analysis undertaken on the calibrators showed 
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Figure 5. Julian-day Intensity Monitoring (JIM) plot of PKS B1519-273 and PKSB1622-253 over the period 2003 - 2005. Large changes 
in the mean flux density, due to intrinsic changes in the intensity of the sources are apparent on timescales of months and years in both 
sources, while the calibrator 3C227 shows no such variations. The larger scatter apparent in the data for PKS B1519-273 and PKS B1622- 
253 compared to the calibrator is primarily due to short timescale ISS-induced variability which is compressed by the time-axis on this 
plot. 
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Figure 8. Changes in the characteristic timescale as a function of day of year for PKSB1622-253 (top) and PKSB1519-273 (bottom). 
The blue circle data points are from the first year of monitoring and the red cross data points are from second year. All-parameter best 
fits shown as solid black lines. Best fits for R=5 shown as dash blue lines. 



weak peaks in the PSD around periods of one day corre- 
sponding to diurnal thermal variations (and also weaker 
harmonics), but these contain much less power than the 
timescales measured for PKSB1519-273 and PKSB1622- 
253 during the same epoch. 

3.1.1 Estimating the error in timescale determination 

Prior to interpreting the apparent changes observed in the 
timescale of variations in PKSB1519-273 and PKSB1622- 
253 it is important to establish the uncertainty in our mea- 
surements. The horizontal error bars in Figure [8] are just 
the length of the observing period. The vertical (timescale) 
error bars were determined empirically from the distribu- 
tion of scintle periods, normalised by the mean period for 
each observing period in order to remove the annual cycle 
effect. This approach assumes quasi-periodicity, as discussed 
above. It also assumes that the statistics of the normalised 
scintle distribution are constant over the whole data set (we 
will show in Section[4]that this is a reasonable assumption). 

Consider a given observing session (typically 10-15 
days) which observes N adjacent scintles with periods Ti 
to Tjv, and a mean scintle period, T perio d, calculated from 
these N values. A period is measured as the scintle peak- 
to-peak or trough-to-trough time. The T perio d estimate is 
associated with the middle of the observing period, with an 
uncertainty in time equal to the observing period. The de- 
viation of scintle period T from T perio d is 

rj-,i -M ^period 

-L dev Ty* 

-t period 

This calculation accounts for the fact that Tperiod may 
change between observing sessions either because of source 
evolution in structure and/or size (which leads to changes in 
the characteristic timescale T c h a r), or because of the annual 
cycle (or a combination of the two) . Figure [9] shows distri- 
butions of Tdev values for adjacent scintles in PKSB1622- 
253 and PKS B1519-273. The distribution observed for both 
sources is very similar and the combined distribution is well 
represented by a normal distribution with a mean of 13.1 
percent and a standard deviation of 6.7 percent. In other 
words, half of all Tdev values calculated from two adjacent 
scintles are less than 13.1 percent, so if an observing period 
has a T c h ar value calculated from two scintles, an error bar of 
0.13 x T c har corresponds to a confidence interval of approx- 
imately 50 percent that the true Tchar value lies within the 
error bars. Assigning a ± 20 percent error bar to the Thai- 
values calculated from two adjacent scintles corresponds to 
a 2-cr (95%) confidence that the true T cria r value lies within 
the error bars. The corresponding error bars for T cria r values 
calculated from 1, 3 and 4 scintles are given by Gaussian 
error theory as 28, 16 and 14 percent respectively. 

The error bars for the timescale axis for each point 
were calculated using the method outlined above, except 
for points where the observed Tdev for an observing period 
exceeded the expected error estimate, in which case the mea- 
sured Tdev was used. 

3.1.2 Fitting an annual cycle model to the data 

The annual cycle in the timescale observed in ISS-induced 
variability can be used to study the scattering screen that 
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Figure 9. The top plot shows the cumulative frequency distri- 
butions of Tdev for each source. The middle plot shows the cu- 
mulative frequency distribution of the combined Tdev values, and 
the bottom plot shows the equivalent frequency distribution. The 
smooth lines in the middle and bottom plots show Gaussian dis- 
tributions of data with the same mean and standard deviation. 



produces it. The relative transverse motion of the Earth 
about its orbit as viewed from the direction of the source 
is fixed for any particular source, however, there are five 
free parameters when fitti ng the observed annual cycle 
ijjauncev fc Macauartll200ll ) : 

V a , Vg: The screen velocity measured in the right ascension 
and declination coordinate system (in km s _1 ) compared to 
the Local Standard of Rest (LSR). 

so: The characteristic scintle spatial scale, is the coherence 
scale (in km) for the turbulence in the screen. The charac- 
teristic scintillation timescale seen by the observer Thar is 
then 



V trans 

where V trans is the relative transverse velocity of the screen 
and observer. 

R, f3: Represents the anisotropy in the spatial pattern of 
the scintles. The scintles are assumed to be elliptical in shape 
with an axial ratio R and the major axis of the ellipse in- 
clined at an angle (3 with respect to the screen velocity vec- 
tor. 

We have undertaken a non-linear least squares fit of 
the annual cycle model to the data for PKS B1519-273 and 
PKSB1622-253. The solid lines in Figure |5] are the optimal 
fits when all parameters are allowed to vary without restric- 
tion. The values of these fits are summarised in Table[T] The 
line of best-fit passes within the error bars of most of the 
measured data points (the error bars showing the 95 percent 
confidence interval) and the reduced \ 2 statistic for the two 
fits are 0.8 and 2.1 for PKSB1519-273 and PKSB1622-253 
respectively. Obtaining a \ 2 value close to 1 for both source 
demonstrates that the error-bar estimation is reasonable. 
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The error estimates in Table Q] were obtained by inves- 
tigating the variation of \ 2 while 2 parameters were simul- 
taneously varied. This approach was taken as the effects of 
the different parameters on the annual cycle fit are not inde- 
pendent, indeed we found the fitting to be quite insensitive 
to the interplay between the anisotropy and LSR velocity 
parameters. The 1-a error estimates given in Table [1] were 
obtained from the width of the \ 2 = 2.3 error ellipse along 
the relevant axis, as is appropriate for systems with two 
degrees of freedom. The interdependence of the different pa- 
rameters in the model means that the annual cycle fitting is 
not a particularly effective way of estimating the anisotropy 
of the turbulent screen. This is also demonstrated in Fig- 
ure [8] where the black line shows the best fit when all pa- 
rameters are allowed to vary and the blue-dashed line shows 
the fit when the anisotropy is restricted to be 7? ^ 5 (the 
best fit for which occurs when R — 5). For both sources the 
best-fit models have an anisotropy ratio more than three 
times greater than this, however, the annual cycles in the 
two cases are not dramatically different. The anisotropy can 
be more directly and accurately estim ated from time-delay 
observations (e.g. iBignall et al.ll2006l) . however, as outlined 
above these are not practical for the more slowly varying 
IDV sources we are dealing with here. As the scintillation 
properties vary with wavelength in a fashion which is both 
significantly different from the differences expected due to 
changes in the source, and predictable, it is clear that simul- 
taneous observations covering multiple frequencies would be 
greatly beneficial. Observations of other sources have also 
seen more extreme variability in the polarized component 
of the flux density, which is interpreted as indicating the 
polarized emission arises from a s maller region than m uch 
of the total intensity emission (e.g. iRickett et alll2002T ). So 
observations at multiple frequencies and preferably also po- 
larization information can greatly aid in disentangle source 
and screen properties and break some of the degeneracies 
present in annual cycle fits at a single frequency. 



4 DISCUSSION 
4.1 PKSB1519-273 

The BL Lac object PKSB1519-273 is listed in the first 
Parkes catalogue, with fl ux densities of 2.0 Jy at 2 . 7 GH z 
and 2.3 Jy at 5 GHz (|Bolton. Shimmins fc Wail 1 1975ft . 
ft is a point source to most instrument s at most wave- 
lengt hs, although both the VLBA at 2- cm ([Lister fc Homanl 
l2005ft and the VSOP 5 GHz survey (|Dodson et all 12008ft 
show it to be slightly extended. There are no confusing 
radio sources nearby, and it is a weak sour ce of soft X- 
rays (|Urrv et all 1 1996ft and possibly of 7-rays (|Fichtel et all 
1 19941 ). Optically, PKS B1519-273 is identified as a B-filter 
magnitude 17.7 star-like obj ect with a featureless spectrum 
|Veron-Cettv fc Verorj|200a ). Visual observations have not 
identified a host galaxy, because it is obscured by nearby 
and partly overlapping objects. PKS B1519-273 was in- 
cluded in high resolution I-band studies of 24 BL Lac ob- 
jects carried out using several telescopes between late 1999 
and early 2001, and found to lie at a redshift of z = 1.294 
based on 4,000 - 10,000 A low-resolution spectral observa- 
tions (|Heidt et al.ll2004ft . A similar redshift was obtained by 



ISbarufatti etafl <|2005l ). In both cases these measurements 
are based on a single emission line in the spectrum identified 
as Mgll and a r e very different from the value of 0.07 given by 
IZensus et~ al. (2002), although it is not clear how or where 
that measurement was obtained. 

PKS B1519-273 is well known to display pronounced 
intrinsic radio variability on time scales of mon ths. The 
ATCA blazar monitoring program of lBignall l|2003ft included 
PKS B1519-273, and eleven sets of flux density observations 
at 8.6 GHz showed it to vary between a mean flux density 
of 1.8 and 2.2 Jy over a 12 month period (starting in Febru- 
ary 2001). Radio variability on time scales of days was dis- 
covered in PKSB1519 -273 at 4.8 and 8.6 GHz duri ng an 
ATCA survey in 1994 (|Kedziora-Chudczer et~al1l200lh . and 
was also obs erved in the 1997-2001 ATCA blazar monitor- 
ing program |Bignal]||2003r i. Multi-frequency ATCA observa- 
tions of PKS B1519-273 were made for five consecutive days 
starti ng on 10 June 1996, an d again starting on 9 September 
1998 (|Macquart et alJl200Cft . 

Multi-frequency ATCA observations of PKS B1519-273 
have found that the transition frequency bet ween strong and 
weak scattering is in the vicinity of 4.8 GHz l|Maccmart et all 
l2000ft . so the COSMIC observations at 6.7 GHz should be in 
the weak scattering regime . Wea k scattering is a broadband 
effect, and lMacquart et alj l|2000ft found that the light curves 
observed at 8.6 and 4.8 GHz were well correlated. There was 
also some correlation in their observations of light curves 
between 4.8 and 2.5 GHz, but there was clearly a significant 
change in the nature of the concurrent scintillation between 
2.5 and 1.4 GHz. 

The annual c ycle i n PKS B1519-273 confirms the work 
of ljauncev et all ( 2003T ) . who concluded that this source dis- 
played an annual cycle in its variability time scale based on 
an examination of ATCA data over ten epochs, including the 
five days in 1998 data; and seven epochs in 2001, each with 
2-3 days of data. Most of these observing periods were too 
short to record more than a single scintle, and the charac- 
teristic variability of each epoch was thus measured as To. 5, 
calculated from the DACF. Three of the ten data points span 
days 120 to 150 of the year, with to. 5 values of ~3 hours (8.6 
GHz) and ~5 hours (4.8 GHz). For quasi-periodic scintles, 
Tchar ~ 6 x to. 5 so at the Ceduna observing frequency of 
6.7 GHz the expected value of T c har is approximately 6 x 
4 hours = 24 hours during this time of year. This agrees 
well with the values observed by Ceduna, as can be see in 
Figure [8] 

4.2 PKSB1622-253 

PKS B1622-253 is listed in the first Parkes catalogue, 
with flux densities of 2. 3 Jy and 2.0 Jy at 2.7 GHz 
and 5 GHz respectively l|Bolton et all I1975T ). It is opti- 
cally very faint and lies behind the Ophiuchus Cloud, 
a dense interstella r gas cloud south of p Ophiuchus 
jHunter et al.lll994ft. PKSB1622-253 has a redshift of z = 
0.786 ( di Serego-Alighieri et al.|[l994ft and an unresolved op- 
tical c ounterpart of B - filter magnitude 21.9 (Schlegel et al, 
1998). iPunslev et~ai1 (120051 ) present VLA radio images of 
PKS B1622-253 made in June 2002, at 3.6 and 6 cm. The 
images show extended structure consistent with the core-jet 
morphology of the central engine model. The main features 
are a powerful radio core, FR II radio lobes consisting of a 
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Table 1. Annual cycle model best-fit values. 



Source 



R 



13 (rad) 



V a (km s" 1 ) V 5 (km s" 1 ) s (10 6 km) 



PKSB1622-253 
PKSB1519-273 



18.0 ±4.0 

15.1 ±3.0 



0.17 ±0.02 
0.07 ±0.02 



-37.5 ±3.0 
-52.0 ±7.0 



8.4 ±0.5 
12.1 ± 1.0 



9.2 ±0.9 
7.1 ± 0.4 



diffuse jet-like eastern extension, and a western lobe that is 
also part of a jet. 7-ray emission from PKS B1622-253 was 
detected by the Energetic Gamma-Ray Experiment Tele- 
scope ( EGRET) on board the Compton gamma ray obser- 
vatory dNolan et al.| [ l99rj ). From the EGRET observations, 
iHunter et al l l| 19941) tentatively determined the 7-ray spec- 
trum of PKSB1622-253 to be described by a power law 
with spectral index -1.9 ± 0.3. However the 7-ray flux is 
PKS B1622-253 is variable, and changed by a factor of ~2 
over the 1991-93 EGRET observ a tion p eriod . 

iTingav. Murphy fc Edwards! (|l998h noted that while 
the Parkes 90 catalogue records PKS B1622-253 as having 
a flux density of 2.02 Jy at 5 GHz, the subsequent Parkes- 
MIT-NRAO survey records its flux density as 3.5 Jy at 4.85 
GHz. Their ATCA observations observed its flux density at 
8.6 GHz to decrease from 2.5 Jy to 0.9 Jy in the 10 month 
period prior to July 1996. Similarly , 22 GHz observations 
with the VLBA ijjorstad et al.ll200ll ) found the flux density 
of the core to vary by a factor of five over the course of a 
year. Radio variability on time scales of days was discovered 
in PKS B1622-253 by the ATCA blaz ar monitorin g pro- 
gram from July 1997 to January 2002 dBignalll I200II ) . The 
source was excluded from the earlier 1994 ATCA survey 
because it failed to meet a compactness selection criterion. 
The source also display s circular radio polarisation variabil- 
ity (|Tingav et al.ll2003l ). but this is too weak to be observed 
by the Ceduna telescope. 

PKS B1622-253 is also assumed to be in the weak scat- 
tering regime at 6.7 GHz. This has not yet been conclu- 
sively determined through multi-frequency observations, but 
the scintillation characteristics of this source are similar to 
PKS B1519-273, and the flux density variability time scales 
of the two blazars follow similar annual cycles. The annual 
cycle in PKS B 1622-253 has not been previously reported 
in the literature. The annual cycle we observe is clearly 
very similar to the annual cycle in PKS B1519-273, although 
the characteristic timescale in this source is typically signif- 
icantly longer. 



4.3 Anisotropy inferred from annual cycle 
measurements 

For both sources the optimum fit to an annual cycle 
implies highly anisotropic scintles and large LSR veloc- 
ity offsets, particularly V a . This is unsurprising as sig- 
nificant velocity offsets are expected for scintillation with 
relatively long variability time scales (days rather than 
hours), which are likely to be associated with more dis- 
tant scattering screens that typically lie hundreds of 
parsecs from Earth, and thus are in motion with re- 
spect to the LSR. Anisotropic scintles have been associ- 
ated with other blazar observations, notably PKS B 0405- 
385 (iRickett. Kedziora-Chudczer fc Jauncevl l2002ri and 
J1819±3845 (jMacauart fc de Bruvnl 120061 ). Many authors 



have argued that, for scintillation due to scattering by the 
ISM, electron density fluctuations are expected to be elon- 
gated along the local direction of the magnetic field (e.g. 
iGoldreich fc Sridhaij 1 19951 ; IChandran fc Backed |2002| ). Al- 
ternatively, anisotropy may be due to an anisotropic source 
structure, however, further observations at multiple radio 
frequencies are required to disentangle the source and screen 
properties. 

The annual cycle effect can be removed from a flux 
density time series by using the best-fit models for annual 
cycles in Tchar to correct the data for departures from a 
specified T^ar value. This enables examination of the over- 
all pattern of scintillation modulations throughout the 2003 
- 05 observing campaign, without the effects of ISM scat- 
tering material motion and anisotropy. This parall e ls the 
exercise carried out by iDennett-Thorpe fc de Bruvnl (120031 ) 
for J1819±3845 data spanning two years of observations. 
Figure [TOl shows all the PKSB1622-253 and PKSB1519- 
273 data for 2003 - 05 with the annual cycle effect re- 
moved by adjusting the data variability to Tchar values of 
4.5 and 3.0 days respectively. The reference Tchar values are 
close to the mean T c har values for the annual cycle, as can 
be seen from Figure [8] Figure | 10l is very similar to Fig- 
ure 12 in iDennett- Thorpe fc de Bruvnl l|2003l ). As was found 
for J1819±3845, the variability of both PKS B1622-253 and 
PKS B1519-273 is, by eye, similar throughout the observing 
campaign. In other words, just as for J1819±3845, the scin- 
tle statistics did not change significantly over the two years 
of observations for PKSB1622-253 and PKS B1519-273. 



5 CONCLUSIONS 

We have detected an annual cycle in the character- 
istic timescale of the variability in the radio emis- 
sion of PKS B1622-253 and confirm the annual cycle in 
PKS B1519-273 determined on the basis of a much smaller 
and more sparsely sampled dataset. The characteristic 
timescale of the variability in these sources ranges from 
around one day (PKS B1519-273 at the fastest time of the 
year) to around 10 days (PKS B1622-253 at the slowest 
time of the year). These findings demonstrate the unique 
range of parameter space being investigated by the COSMIC 
project. Such long characteristic timescales cannot be prac- 
tically monitored using premier interferometer systems such 
as the VLA or ATCA because these instruments are in too 
high demand to dedicate long periods of time to flux density 
monitoring. The best-fit annual cycle for both these sources 
implies scattering screens with significant transverse veloci- 
ties with respect to the LSR and a high degree of anisotropy. 
Through the COSMIC program, we will continue to mon- 
itor these and other strong southern IDV sources to look 
for signs of changes in the source structure and scattering 
screen. 
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Figure 10. PKSB1622-253 (top) and PKSB1519-273 (bottom) data adjusted by the time scales calculated from the best fit annual 
cycle models to remove the annual cycle effect, with the modified time series adjusted to "constant" variability timescales of 4.5 days 
and 3.0 days respectively. 
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